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ABSTRACT 



In network and active region plages, the magnetic field is concentrated into structures often described as flux tubes (FTs) and sheets 
(FSs). Three-dimensional (3-D) radiative transfer is important for energy transport in these concentrations. It is also expected to be 
important for diagnostic purposes but has rarely been applied for that purpose. Using true 3-D, non-local thermodynamic-equilibrium 
(non-LTE or NLTE) radiative transfer (RT) in FT and FS models, we compute iron line profiles commonly used to diagnose the 
Sun's magnetic field by using and comparing the results with those obtained from LTE or one-dimensional (1-D) NLTE calculations. 
Employing a multilevel iron atom, we study the influence of several basic parameters such as either FS or FT Wilson depression, 
wall thickness, radius/width, thermal stratification or magnetic field strength on Stokes / and the polarized Stokes parameters in the 
thin-tube approximation. The use of different levels of approximations of RT (3-D NLTE, 1-D NLTE, LTE) may lead to considerable 
differences in profile shapes, intensity contrasts, equivalent widths, and the determination of magnetic field strengths. In particular, 
LTE, which often provides a good approach in planar 1-D atmospheres, is a poor approximation in our flux sheet model for some of 
the most important diagnostic Fei lines (524.7 nm, 525.0 nm, 630.1 nm, and 630.2 nm). The observed effects depend on parameters 
such as the height of line formation, field strength, and internal temperature stratification. Differences between the profile shapes may 
lead to errors in the determination of magnetic fields on the order of 10% to 20%, while errors in the determined temperature can reach 
300 - 400K. The empirical FT models NET and PLA turn out to minimize the effects of 3-D RT, so that results obtained with these 
models by applying LTE may also remain valid for 3-D NLTE calculations. Finally, horizontal RT is found to only insignificantly 
smear out structures such as the optically thick walls of flux tubes and sheets, allowing features as narrow as 10 km to remain visible. 

Key words. Line: formation - Magnetic fields - Polarization - Radiative transfer - Sun: atmosphere - Sun: photosphere 



1. Introduction 

The magnetic fields in both the solar network and active region 
plages is concentrated into elements that are generally described 
by magnetic fl ux tubes and sheets (IStenflolll973l ISpruillll976t 
ISolankil[l993l) . These magnetic elements form the footpoints of 
either coronal funnels and loops and chromospheric canopies. 
As such, they channel the energy tra nsported from the solar 
interi or to the upper solar atmosphere dNarain & Ulmschneiderf 
1 19961) and play an impo rtant role in modulat ing the Sun's total 
and spectral irradiance dDomingo et al.ll2009l) . 

Magnetic elements are thought to be formed by the magneto- 
convective processes of flux expulsion and convective collapse. 
Flux expulsion sepa rat es the magn e tic flux from the convec- 
tion dParkerfll963allbL IWeissI 119661: iGallowav & Weiss! Il98lt 
Hurlburt & Toomre 1988]), while c onvective colla p se bring s field 
strengths up to kG levels (e.g. Parker] 
Grossmann-Doerth et al 1998l lLites etalT 
2010). The magnetic elements have a range of sizes down to 
as small a s 100 km, close t o the highest currently achievable 
resolution dLagg et al.l l2010). and probably lower than that too 
(Riethmiiller et al., in prep.). 

Radiation plays a dominant role in determining the proper- 
ties of magnetic elements, both during the formation — i.e. dur- 
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ing the convective collapse, where radi ation leads to the cooli ng 
of the gas in the magnetic features, cf. lVenka takrishnanl d 19861) — 
and during the normal lifetime of a magnetic element. Since the 
kG magnetic field suppresses the convection within the mag- 
netic elements, their photospheric layers are heated mainly by 
radiation. Owing to the evacuation produced by the excess mag- 
netic pressure within the magnetic elements, radiation flows in 
throug h the walls a nd heats the interiors of the magnetic ele- 
ments dSpruitl fl976). Multidimensional radiative transfer (RT) 
plays a central role in setting this temperature and is a standard 
part of modern MH D simulations of magnetoconvection in the 
Sun's su rface layers dNordlundll 19821 11983L 11985b iDeinzer et ail 
1984a b; Vode r et al.ll2004l 120051) . 



Diagnostic RT, however, is generally carried out in 1-D or 
at most in so-called 1.5-D, i.e. along many parallel rays pass- 
ing through a multidimensional model (e.g. iBtinte et al.l fl993: 
ISolanki| [l989). This difference in the treatment of the energetics 
and the diagnostics is a basic inconsistency that is worth it to be 
investigated. We henceforth also use the term 1-D for 1.5-D RT. 
Although there have been a n umber of 3-D RT studies of convec- 
tive features dCarlssonll2009l and references therein), the multi- 
dimensional RT diagnostics of magnetic elements are to our 
knowledge restricted to only a few pape rs dStenh olm & S tenflol 
ll977lll978l:rBruls & von der Liihe 2001). 
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Owing to the very restricted computer capacities at the time, 
Stenholm and Stenflo in their pioneering papers had to make 
some compromises, restricting themselves to a strictly cylin- 
drical flux tube and a two-level atom. In the present work, 
we use a more sophisticated (but still simple) flux-tube model 
where its expansion with height is computed co nsistently using 
the thin-tube (t hin flux-sheet) approxi mation dLeightonl 1963: 
lParkerl!974allrxlSpruit & Roberts! 19831 and references therein), 
a multi-level atom, and true 3-D NLTE RT. A comparison with 
3-D radiation MHD simulations has shown that the thin-tube 
(or thin flux-sheet) approximation provides a reasonable de- 
scription of the magnetic flux concentrations in the simulations 
(lYelles Chaouche et al.ll2009l) . 

The thin FS geometry was studied using 2-D RT by 
iBruls & von der Liihd d200l | ) and part of our work resem- 
bles theirs. However, while IBruls & von der Liihd d200lh aimed 
mainly at establishing that features narrower than the horizon- 
tal photon mean free path can still be resolved, the aim of the 
present paper is broader. We aim to quantify the influence of 
multi-dimensional NLTE diagnostic RT relative to LTE and 1-D 
NLTE RT. To this end, we carry out a parameter study probing 
the influence of magnetic field strength, flux-tube radius, flux- 
tube atmospheric model,and UV opacity on the Stokes profiles 
of four commonly used spectral lines. Although MHD simula- 
tions have reached a high level of realism, we prefer to be able 
to change parameters at will in a simpler model. This allows us 
to get a clearer understanding of the effects of multi-dimensional 
RT in comparison with 1-D NLTE and 1-D LTE. 3-D diagnos- 
tic RT in hydrodynamic and MHD simulations will be treated in 
forthcoming papers. 

The structure of our work is as follows: In section|2l we pro- 
vide a summary of our model ingredients. In section [3] we de- 
scribe the results of our parameter study in detail. In Section |H 
we discuss the influence of the parameter variations on widely 
used quantities, such as contrasts and equivalent widths, fo- 
cussing on the differences between three calculation methods 
(LTE, 1-D NLTE, and 3-D NLTE). Our conclusions are pre- 
sented in section [5] 



2. Model ingredients 

2.1. Radiative transfer calculations 

All our calculations were done with the RT code RH 
adapted to our ne eds. This versatile code was developed by 
Uitenbroekl d2 000) and is based on the iteration scheme of 



Table 1. /-values of our four selected lines 



Rybicki & Hummer l (119911 119921, |1994l) . It uses an improved 



version (ISocas- Navarro et al. 2000|) of th e sho rt characteristics 
(SC) m ethod of lOlson & KunaszM 19871) and iKunasz & Auerl 
(1988) for Stokes I and lRees et al.ld 19891) for the full Stokes vec- 
tor. The code is able to do LTE or NLTE calculations in 1 -D, 2-D, 
and 3-D. If not otherwise stated, the calculations pres ented here, 
were always done with the 3-D version of the code dUitenbroekl 
2006). In addition to several adaptions for convenience or to en- 
hance the performance, we made a principal change to the for- 
mal solver that has a considerable influence on our results. In the 
appendix, we summarize our improved method for the formal 
solution used for both Stokes / and full Stokes vector calcula- 
tions, as well as some pitfalls we encountered with both the old 
and new methods. 
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2.2. Atomic model and iron abundance 

We used an Fei atom including 22 Fei levels and the ground 
state of Fe n, which allowed us to perform a realistic NLTE cal- 
culation of selected multiplets containing transitions of interest. 
The model, kindly provide d by J. H. M. J. Bruls, employs en- 
ergy level data compile d bv lCorliss & S ugar ( 1982) and line data 
from lFuhr etall d 1988b. It was en hanced with oscillator strength 
data from lTheveninl(ll989[ll990l) . 

We selected the well-studied and often observed 524.7 nm, 
525.0 nm, 630.1 nm, and 630.2 nm lines for this w ork, which are 
often analyzed eithe r individually o r as line pairs (|Stenflol 
Elm ore et all f l992: Ts uneta et all 120081: ISolanki et alT 
Martinez Pillet et al.ll201 lb . Altogether, the atom contained 23 
continuum and 33 line transitions with a grand total of some 
1300 wavelength points (we tried to keep this number small to 
reduce the computing resources needed). In this model atom, the 
oscillator strengths were as displayed in Table Q] The /-values 
of the 524.7 nm and 525.0 nm lines were enhanced by roughly 
a fac tor of from two to thr ee against those found in the litera- 
ture dBlackwell et al. 19791). follow ing the correction procedure 
used in Bruls & von derLuhdf2 001). An iron abundance of 7.44 
was employed following lAsplund et al.l d2000h . We note that we 
were mainly interested in differences between the line profiles 
obtained in LTE, 1-D NLTE, and 3-D NLTE under different con- 
ditions and not in the absolute line profiles. In particular, we did 
not compare our line profiles with observations, so that details 
such as the exact elemental abundance or /-values are not of 
central consequence and do not influence our conclus ions. 

T he well-known problem of UV overionization dBruls et al.l 
119921) . which leads to wrong population numbers because of 
missing iron line opacity in the UV, was eliminated by enhanc- 
ing the opacities in the relevant wavel ength range (opacity fudg- 
ing) as proposed by Brul s et aTl dl992l) . The resulting differences 
were calculated and are discussed in Sect. 13.21 



2.3. Model atmospheres and computational setup 

We studied both magnetic flux tubes (displaying axial symme- 
try) and flux sheets (displaying translational symmetry). In both 
cases, we employed the thin-tube approximation to describe 
their internal structure. For most computations, however, we 
used magnetic fields only to determine the Wilson depression 
and the expanding geometry of the flux t ube, to better isolate 
the infl uence of multidimensional RT, as in Stenholm & Stenflol 
(fl977h . In Sect. 13.91 however, we also consider results computed 
by including the formal solution of the Unno-Rachkowski equa- 
tion, i.e. where w e computed Stokes profiles using the so-called 
field-free method dReeslll969l) . Thus, we included the influence 
of the Zeeman effect on the line profiles (splitting, polarization), 
although the results presented here neglect the influence of the 
Zeeman effect on the population numbers. 

The above-mentioned flux sheet could in principle be calcu- 
lated with a 2-D code as long as one considers radiation emerg- 
ing only along rays lying in a plane perpendicular to the plane 
of symmetry. However, to achieve an unbiased comparison with 
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Fig. 1. Temperature distribution in the atmospheric cross- 
section through a flux sheet model. The t c = 1 level (thick line) 
and the nominal boundaries of the flux sheet (thin vertical lines) 
at t c — 1 are indicated. 



flux tubes, we used the 3-D code version. The qualitative effects 
in a 2-D flux sheet (FS) are the same as in a true 3-D flux tube 
(3-D FT), but are achieved at much lower computing costs. We 
note that an equidistant grid (in the xy-direction enforced by the 
code) is not very well-suited to mapping an expanding tube or 
sheet, therefore the number of grid points has to be quite high. 
To achieve reliable results, we used models with 288 x 4 x 256 
grid points in the flux sheet geometry, where the z-axis denotes 
the height in the atmosphere (see also Fig.[T]for a cut across the 
flux sheet). In Sect. 13.71 we present the results for a true 3-D flux 
tube geometry having 144 x 144 x 128 grid points in the carte- 
sian geometry of the RH code. Although only half the number 
of pixels are introduced along the x and z axes when compared 
to the flux sheet case, the total number of grid points is an order 
of magnitude larger. To keep the computing time from growing 
accordingly, the number of rays was reduced (standard set A6 
of RH code instead of 54 vertical directions per octant, see Sect. 

531 , 

To compare our 3-D with our 1-D calculations, we con- 
structed plane-parallel atmospheres by vertically cutting the 3-D 
cube at selected x-positions. The 1-D NLTE RH code was then 
run with the so-constructed 1-D atmospheres as input (we stick 
to the term 1-D NLTE for this method, which is often called 
1.5-D NLTE). In addition, for comparison purposes 1.5-D com- 
putations in LTE were carried out for each flux sheet or tube 
(henceforth simply called LTE computations). 

2.3.1. Standard model 

We varied a range of parameters such as radius, Wilson depres- 
sion, thickness of the walls, filling factor, and flux sheet atmo- 
sphere (see Sect.[3]for details). As a reference, we defined a stan- 
dard model with a fixed parameter set: 

Atmosphere: If not otherwise stated, the inner and outer 
atmosp heres were analogous to those in IStenholm & St enflo 
(Il977h and therefore identical. This choice allowed us to more 
clearly segregate the effects of horizontal RT. We used the 
Harvard-Smithonian reference atmosphere (HSRA) described 



bv iGingerich et al] (Il97ll) . In the field-free part, the model was 
extended downwards by combining it with the convection zone 
model of ISpruitl (|1974|) . Figure Q] shows the geometry and the 
temperature distribution in the xz plane for our standard model. 
The thick solid line follows r c = 1 and outlines the smooth but 
quite steep transition between the flux sheet interior and exte- 
rior. Inside the nominal boundary (indicated by the vertical lines 
as in many subsequent figures), the hot walls of the outer atmo- 
sphere, positioned below the continuum level of the quiet Sun, 
lie close to the local t c = 1 level, so that their presence strongly 
influences the radiation field sensed by an atom within the mag- 
netic feature. A constant microturbulent velocity of 0.8 km/s was 
added to every grid point of the model in direct correspondence 
with the original model. In the following, for simplicity, we refer 
to the outer, field-free atmosphere as either the quiet Sun (QS) 
or QS part. We note that the outer atmosphere ceases to exist at 
the height where the canopy fills the whole width of the atmo- 
sphere except for one single pixel at the edge of the atmosphere 
that was chosen to be pure HSRA without any magnetic field. 

Field strength: The magnetic field strength of our standard 
sheet was B = 1600 G at r? s = 1 and B = G outside, re- 
sulting in a Wilson depression of 191 km. This B value can- 
not be directly compared with measured intrinsic field strengths, 
which depend on the height of formation of the employed spec- 
tral line(s), diagnostic method, etc. The field strength dropped 
with height and consequently the magnetic element expanded 
with height. We note that a flux sheet exhibits a more rapid ex- 
pansion with height than a round flux tube as a consequence of 
flux conservation if one assumes a fixed run of field strength 
with height because B x R — const, for a sheet in contrast to 
BxR 2 = const, for a tube, where R is the radius of the FT or the 
horizontal location of the FS boundary. 

Wall thickness: The sheet boundaries were smoothed in the 
x-direction with a Fermi function which drops from 90 to 10 
percent within 20 km, that is sli ghtly above the m inimal value of 
the wall thickness estimated by Schussler ( 1986). 

Radius, filling factor, and grid points: Our standard model 
had a radius of R = 100 km and a filling factor a = 0.2. The grid 
points were equidistant in each dimension, where dx = 3.47 km, 
dy = 100 km (irrelevant for a 2-D geometry), and dz = 5 km, 
starting at 460 km below t c = 1 of the field-free part of the 
model. Such a low starting point was needed to ensure that the 
model is optically thick at the bottom, even in the presence of a 
strongly evacuated flux sheet. 

2.3.2. Angular discretization 

Multi-dimensional problems basically need higher angular res- 
olution to achieve the same accuracy o f angular quadrature as 
in a 1-D problem (Ivan Noort et a 1. 2002). The very specific ge- 
ometry of FS/FTs compelled us to increase the number of rays 
used in the simulation. We used a grid of 27 vertical (i.e. inclina- 
tion) and 2 horizontal (i.e. azimuthal) angles per octant in the FS 
models, i.e. 216 outgoing rays from each grid point. This small 
number of horizontal angles was adequate in the (2-D) FS ge- 
ometry. Figure [2] illustrates the difference in outgoing intensity 
between our chosen angular resolution and one obtained with the 
often used A6 standard set of the RH co de with only 24 outgo- 
ing rays (cf. lBruls & von der Close to the flux sheet 
boundary, the difference is significant (« 7%). This difference 
has its origin in the very hot walls of the FS in this model and 
the optical thinness of the sheet's interior. The part where the 
walls are really hot is very localized, being restricted to a few 
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Fig. 2. Vertically emergent intensity in the core of the 525.0 
nm line as a function of x-position, counted from the FS center 
(i.e. from the horizontal axis of Fig. Q] cutting across the FS). 
The plotted curves were calculated with two different angle sets, 
where solid lines indicate the set of 27 (inclination) x 2 (az- 
imuth) rays per octant and dashed lines the standard set A6 of the 
RH code with 6 rays per octant. The thin vertical line indicates 
the nominal FS boundary (R = 100 km), which was reached at 
the t c = 1 level in the field-free atmosphere (at the same position 
as the rightmost vertical line in Fig. Q}. 



grid points. Because the strongly evacuated sheet interior is op- 
tically thin, the hot walls can be seen from far away, as long as 
the ray travels only upwards through the sheet interior. The dis- 
cretization of the rays gives rise to different lighthouse-like rays 
through the sheet. To see the hot walls at every point within the 
sheet it was necessary to increase the number of vertical rays. 
Some locations in the interior would not otherwise have been hit 
by a ray originating in the hot wall, which, in turn, may affect the 
local atomic level populations and the corresponding line depths. 
We note, however, that the artificially suppressed NLTE effects 
at these points are only relevant when the optical thickness of 
the interior is not negligible. 

For the 1-D NLTE calculation, we used a standard set of five 
Gaussian angles for all calculations. 



3. Results 

3.1. General features 

We like to compare our results for the LTE, 1-D NLTE, and 3-D 
NLTE RT computations. We first concentrate on the Fei 525.0 
nm line formed in our standard atmosphere (see Sect. 12.3. U . but 
later also consider the other spectral lines and models. Figure [3] 
shows the intensity along the x-axis of the standard flux sheet 
(i.e. cutting across the sheet) for three selected wavelength posi- 
tions, namely in the line core, in the line flank at approximately 
half depth (~ 0.004 nm from the line center), and near the contin- 
uum (~ 0.008 nm from the line center). The nominal flux sheet 
boundary at r c = 1 is indicated by vertical lines. The center of 
the flux sheet is located at x = 0. Only half of the cut through 
the sheet is plotted since the other half is simply a mirror image 
owing to the symmetry of the model. 
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Fig. 3. Intensities calculated with the different methods, indi- 
cated in the top panel, as a function of x position at three differ- 
ent wavelengths in the 525.0 nm line: a) line core; b) half depth 
(~ 0.004 nm); c) near continuum (~ 0.008 nm). The nominal FS 
boundary (R = 100 km at t c = 1) is indicated by the vertical 
line. 



In the line core and flank, both 1-D and 3-D NLTE compu- 
tations sh ow similar, strong variations along x (Fig. |3j, as re- 
ported by Bruls & von der Luhek 200R with the rest intensity 
reaching a peak value of two to four times that found in the QS 
far away from the FS. We note that quantitative differences be- 
tween the results of iBruls & von der Liihd ([2001) and ourselves 
may be caused by the quite different atomic and atmospheric 
parametrizations, the most important probably being the Wilson 
depression (magnetic field strength). In the line wing, the vari- 
ation is slightly higher than in the line core. The LTE intensity, 
however, shows only a minor variation along the x-axis in both 
the line wing and core. Far away from the FS, all three meth- 
ods of computation give almost identical results, only the LTE 
rest intensity remaining slightly below that of the NLTE meth- 
ods. In the continuum, all three methods are identical and show 
a s harp rise in the rest intens ity close to the boundary of the FS 
(cf. lCaccin & Severinol[T979h . 

Finally, a large difference between the 1-D and 3-D NLTE 
calculations is found inside the FS, with the 3-D profile being 
weakened (i.e. having a higher rest intensity than the profiles 
computed with the other methods). 

The upper panels of Fig. |4]show the corresponding spectra 
for five selected locations indicated above the panels. In all plot- 
ted cases, the LTE line profile saturates, maintaining nearly the 
same rest intensity (as seen in Fig. [3^, while the NLTE profiles 
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Fig. 4. Intensity profiles of the 525.0 nm line at five selected x positions as indicated at the top of each panel (see also Fig.[3]l. Solid 
thick line: 3-D NLTE; dashed line: LTE; plus signs: 1-D NLTE. Upper panels: with opacity enhancement according to Brul s et al.l 
( 1992); lower panels: no additional UV opacity. 



become much weaker at some locations. We note that the rel- 
ative differences in the line depth (see Fig. HJ are considerably 
smaller than those in rest intensity (Fig.O owing to the generally 
low line core intensity. 

We now attempt to understand the intensity variations plot- 
ted as a function of x in Fig. [3] Far away from the flux sheet, 
the line is formed in the QS part of the atmosphere just below 
the optically thin canopy, which lies at a height of 400 km. The 
difference between the LTE and NLTE profiles at that location 
is basically the same as one gets when using the HSRA on its 
own, without any FS. As we move closer to the sheet, the height 
of the canopy's base decreases continually. Unit optical depth in 
the line is reached at lower heights (below the canopy and, for a 
sufficiently low canopy base, at its boundary) owing to the lower 
density and hence opacity of the canopy, which remains optically 
thin to the line right up to positions slightly inside the nominal 
FS boundary. The line core is therefore formed in the hotter lay- 
ers of the external atmosphere and consequently has a higher line 
core intensity, which culminates just inside the sheet's nominal 
boundary. 

The 3-D NLTE contribution functions in Fig.|5^ illustrate the 
decreasing height of the line-core formation region as we move 
towards the sheet. We note the peculiar shape of the contribution 
functions far away from the sheet. In contrast to the typical case 
in which the contribution function first sharply increases when 
moving up in the atmosphere, then slowly decreases above the 
height of maximum contribution, these contribution functions 
show a steep edge at the higher end, which corresponds to the 
lower boundary of the canopy. Examples are found between the 
two highlighted contribution functions (solid lines) at x = 465 
and x = 150 km. 

The position of this edge moves down as we move closer to 
the sheet in complete accordance with the geometry of our at- 
mosphere (see also Fig. Q]). The intensity maximum of I core for 
the NLTE-cases lies slightly inside the flux sheet, where we — 
looking down vertically — directly see the hot walls below the 
t c = 1 level of the QS part. All three contribution functions (at 
the line core, at half-depth, and in the continuum, i.e. Fig.|5ji,[5]3, 
and[5j;, respectively). At that position, the contribution functions 



are peaked almost exactly at the same height, confirming the hy- 
pothesis that we see light formed below the height at which the 
continuum is formed far away from the FS. We note that the 
spectral line does not completely disappear owing to a small dif- 
ference between the line and continuum formation levels and the 
very strong temperature gradient at the hot walls. 

Why do we not see the same effect in the line formed under 
LTE? As the population numbers of Fe i are generally higher in 
LTE (owing to missing ionization by UV radiation), the line is 
slightly stronger than in the NLTE case. This leads to the situa- 
tion where the canopy starts to become optically thick in the line 
core even far from the FS. Closer to the FS, the optical thick- 
ness of the canopy increases as the height of its base decreases. 
Because the temperature in either the canopy or the FS at heights 
of around 300 - 400 km is only slightly cooler than in the ex- 
ternal atmosphere at the same geometrical height, the line core 
intensity decreases only slowly when we move towards the sheet 
(i.e. the line becomes slightly stronger). 

Moving away from the line core, the opacity decreases, so 
that in the LTE continuum the intensity shoots up, as in the NLTE 
line cores (Fig. [3]: ). The three identical continuum curves in Fig. 
[3J; indicate that the continuum is as expected formed in LTE. In 
the line flanks, the LTE opacity takes on an intermediate value 
so that the intensity still increases when horizontally approach- 
ing the FS boundary until reaching a gentle maximum at the FS 
boundary (see Fig.[3j> ). 

The main effect of 3-D RT is seen within the sheet. There the 
LTE and 1-D NLTE profiles attain nearly the same depth as they 
have far from the flux sheet, as expected for a model with the 
same atmospheres inside and outside the flux sheet. Inside the 
flux sheet the LTE and 1 -D NLTE profiles are independent of x 
because the local atmosphere does not change (when neglecting 
the wall smoothing). The 3-D NLTE line profile, however, re- 
mains significantly weaker (cf. Fig. [4^) owing to the ionization 
of neutral iron by the UV radiation flowing i n sideways from the 
hot w alls of the flux sheet, as discussed bv lStenholm & Stenflol 
(Il977h . The line-core contribution function becomes very broad, 
extending from the continuum forming level at -200 km up to 
+300 km in the sheet. 
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Fig. 5. Contribution functions computed in 3-D NLTE at differ- 
ent ^-positions for the same three selected wavelengths of the 
525.0 nm line as in Fig. [3] a) line core, b) line flank, c) contin- 
uum. The dotted contribution function is given for every tenth 
grid position corresponding to Ax « 35 km, starting at 10 km 
from the plane of symmetry of the flux sheet. The five contribu- 
tion functions corresponding to the same x-positions as the pro- 
files shown in Fig. |4] are highlighted (solid lines) and indicated 
by their x-position. Note the very different shapes and positions 
of the individual curves. 



In a small region at the FS boundary, the 1 -D NLTE intensity 
exceeds the 3-D case. In the hot walls, both the horizontal RT to- 
wards the cool sheet interior and radiation from the cool interior 
of the flux sheet (which are both not present in the 1-D case) de- 
crease the amount of UV radiation reaching an atom there, thus 
giving rise to the opposite effect with respect to that seen inside 
the flux sheet. 

We conclude this section by restating Sten holm & Stenflol 
(fl977h since they provided a succinct summary of the results 
presented so far: "Even when a line in a plane-parallel atmo- 
sphere is formed in LTE (i.e. LTE is a reasonable approxima- 
tion), there is no guarantee that LTE is still a good choice in 
more complex geometrical configurations." 



3.2. The role of UV opacity 

It is important to judge th e influe nce of enhancing the UV opac- 
ity followi ng Bruls et al. J 19921) to compare with the work of 
IStenholm & Stenflol d 1977b . who used a simple two-level atom 
without any enhanced UV opacity. The role of UV opacity is 
to partially suppress the influence of the UV radiation from the 
deeper, hotter layers and therefore reduce the NLTE effect of U V 
overionization. 

The spectra in the lower panels of Fig.|4]are calculated with- 
out any enhanced UV opacity. Even far away from the sheet, the 
NLTE-profiles differ rather significantly from the LTE profiles 
when no opacity enhancement is introduced. They start to differ 
dramatically when approaching the sheet (at 150 km). The dif- 
ference between the 1-D and the 3-D calculations in the sheet is 
also strongly enhanced due to the strong influence of the hor- 
izontal transport of UV radiation in the 3-D case. The depth 
of the line in the 3-D case is smaller by more than a factor 
of two. These results are mo re comparable to those obtained 
bv lStenholm & Stenflol (1 19771) . Obviously the employment of a 
more adequate UV opacity model reduces not just the strength 
of NLTE effects in general, but in particular the effect of the hor- 
izontal RT, which remains important nontheless inside the mag- 
netic element. 

We note that the influence of the UV opacity enhancement 
on the Fei 630.1 nm line, which is principally more sensitive 
to NLTE effects, in particular horizontal transport, is slightly 
smaller than that on the 525.0 nm line. The 630.1 nm line is 
already in a partially saturated regime and therefore removing 
the UV-opacity enhances the NLTE effects to a lesser extent. 

3.3. Behaviors of different lines 



The discussion in Sect. 13.21 already indicates that the effects of 
NLTE and 3-D RT depend strongly on the formation height, so 
that different lines may exhibit different behaviors. In this sec- 
tion, we compare the profiles of three widely employed lines 
(Fei 525.0 nm, Fei 630.1 nm, and Fei 630.2 nm). The results 
for the 524.7 nm line turn out to be exactly the same as for the 
525.0 nm line, which is unsurprising, given that the two lines are 
nearly identical, save their Lande factors (which do not play a 
role here, since we set B — G for the line-formation computa- 
tions). The 630.1 nm line is expected to be somewhat stronger, 
and the 630.2 nm line to be somewhat weaker than the 525.0 
nm line. In addition, the two 630 nm lines have a considerably 
higher excitation potential, which may also influence our results. 

Figure [6] shows the line core intensities along the x-direction 
as in Fig. [3^, but now for the three lines 525.0 nm, 630.1 nm, 
and 630.2 nm. We repeat the results for the 525.0 nm line from 
Fig- Eh for ease °f comparison. The global features are the same 
for all three lines. Far away from the sheet, the lines are formed 
in a very similar way (although the 630 nm doublet is somewhat 
stronger in NLTE relative to its strength in LTE). Close to the 
flux sheet boundary, there are, however, notable differences. 

A major difference of the 630.2 nm line in comparison to 
the other lines is the generally higher intensity level owing to 
its weaker opacity. The LTE curve confirms this view, as can be 
deduced from the increasing I core level as we move towards the 
sheet: this is because the LTE line is formed (partially) in the hot- 
ter gas below the canopy and not mainly either in the canopy it- 
self or its boundary where the cores of the other lines are formed 
(see also the discussion in Sect. 13.2b . The significantly stronger 
3-D effect seen in this line is quite striking inside the FS com- 
pared to that for the 525.0 nm line. 
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Fig. 6. Line-core intensity of three spectral lines (identified above each frame) along the x-axis of the model. Solid: 3-D NLTE, 
dashed: LTE, plus signs: 1D-NLTE. 
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Fig. 7. Line-core contribution functions as in Fig. [5^ for the 
three lines 525.0 nm (a), 630.1 nm (b), and 630.2 nm (c). 



The 3D-NLTE case of the 630.1 nm line differs quite 
strongly from those of the other lines. Although the intensity 
increases towards the sheet, there is no peak near the flux sheet 
boundary, which remains at almost the same level over the whole 
sheet. The line-core contribution functions drawn for the three 
lines in Fig. [7] indicate that the weaker lines have strong contri- 
butions from the continuum near the boundary (at x = 94 km), in 
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Fig. 8. Line-core intensity of the Fei 525.0 nm line along the 
x-axis for three different atmosphere smoothings (thick curves 
40 km; normal 20 km; thin 10 km). 3-D NLTE solid curves on 
the left; 1-D NLTE solid right; LTE dashed. 



contrast to the stronger 630. 1 nm line. The contribution function 
in the line core is clearly shifted towards higher levels, i.e. we do 
not see directly the hottest parts of the walls. We note that out- 
side the sheet (e.g. at x = 150 km), the 630.1 nm line is formed 
only marginally above the other lines. Inside the sheet, the effect 
of horizontal RT weakens this line, leading to higher intensities 
in the 3D-NLTE case. 

Figure |6j5 indicates that at certain x-positions the LTE in- 
tensity of the 630.1 nm line is slightly higher than in NLTE. 
This may be explained by mainly two different effects (Jo Bruls, 
private communication). The photospheric source functions for 
these lines usually follow the Planck curve for both LTE and 
NLTE. In the NLTE case, owing to the lower population of Fe I, 
the r = 1 level is shifted downwards, into hotter regions, weak- 
ening the line. On the other hand-especially for the stronger 
630. 1 nm line formed higher up-the NLTE source function may 
clearly fall below the Planck function, leading to the opposite 
effect. Which effect wins in the end is determined by the com- 
plex interplay of many transitions and levels in the atom and is 
difficult to forecast. 



3.4. Influence of the wall thickness 

The true thickness of the walls of magnetic elements is not well- 
known. The minimum thickness based on the expe cted solar 
magnetic resistivity, diffusivity etc. was estimated by ISchtissleil 
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(1986) to be on the order of 3 - 10 km. In addition, radiation 
is expected to smooth the transition from the surroundings to 
the FT/FS and convection is expected to make the wall irregu- 
lar, although the profile shape of the g — 3 line at 1 .56^/m sug- 
gests that the magnetic field drops rather suddenly at the edge 
of the magnetic element dZaver et al.ll 1 989h . In our calculations, 
we used smoothings with three different Fermi functions, each 
defined by the distance over which one atmospheric component 
changes from the 90% to the 10% level. In total, thicknesses of 
40 km, 20 km (standard model), and 10 km were tested. We ex- 
pected the wall thickness to have a strong influence on the lines 
where the continuum radiation of the hot wall could be seen di- 
rectly, as in the case of the 525.0 nm or the 630.2 nm lines. By 
changing the distance over which the atmosphere changes from 
the cold interior to the hot QS outside the sheet, we directly 
changed the temperature at the locations where the continuum 
radiation forms. We also expected that the effect of horizontal 
RT would decrease when we expanded the thickness of the flux 
sheet boundary, since a smoother transition brings the FS ever 
closer to the plane-parallel, i.e. the 1-D case. The resulting in- 
tensity is determined by a complex interplay of many factors 
(see also iBruls & von der Luhell200ll) . 

Figure[8]shows the LTE, 3-D NLTE, and the 1-D NLTE line- 
core intensities of the 525.0 nm line along the x-axis for the 
three different wall thicknesses. As expected, the strongest influ- 
ence is seen at the hot wall itself. While the LTE curves almost 
coincide, the 1-D and 3-D NLTE curves show a systematic in- 
crease in the intensity maximum with decreasing wall thickness. 
Interestingly, the change of the 3-D intensities in the inner part 
of the sheet is, however, much smaller. The largest effect is dis- 
played by the 524.7/525.0 nm lines. The other lines (not shown), 
especially the 630.1 nm line, are less sensitive to the wall thick- 
ness. 
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Fig. 9. Intensity in the line core of the three lines (525.0 nm 
upper, 630.1 nm middle, and 630.2 nm lower panel) at selected 
x positions (diamonds maximum intensity near boundary of FS, 
plus signs FS center, asterisks QS far away from FS) for models 
with different half- widths R of the FS. At R = 100 km, we have 
added the corresponding LTE values, which do not vary with FS 
width. 



3.5. Variation in the radius 

We next considered FS models with different widths. To change 
the width, the flux sheets were stretched in the x direction. 
Therefore, the thickness of the walls scaled with the width of 
the FS. All other quantities remained the same. 

Core intensities at three locations along the x-axis are plot- 
ted in Fig. |9]versus the FS half-width of three spectral lines. The 
selected locations are the sheet center (plus signs), the position 
of the intensity maximum near the wall of the FS (diamonds), 
and the QS far away from the sheet (asterisks). The latter does 
not show much variation with radius as a consequence of its po- 
sition far away from the sheet. For r > 400 km, the intensity 
at the FS center almost corresponds to the 1-D NLTE case. As 
expected, these broad flux sheets are optically thick in the hori- 
zontal directions, so that the hot walls are no longer visible at the 
center of the FS, hence do not influence the radiation field there. 

The intensity maximum near the boundary (diamonds) dis- 
plays a similar behavior. This can be understood by bearing in 
mind that the wall thickness increases with FS width, so that for 
wider tubes the walls are not so hot (see Sect. 13.41 1. However, 
horizontal transport plays a role as well, since an observer lo- 
cated close to one wall also sees radiation from the opposite wall. 
With increasing FS half-width, the strength of this irradiation de- 
creases. For the intensity maximum at the flux sheet boundary, 
the effect of horizontal RT seems to saturate at R « 200 km. This 
is consistent with the behavior seen at the center of the flux sheet 
where the effect saturates at double the distance, i.e. at R « 400 
km. Therefore, the opposite wall (at 2R) is invisible at the bound- 
ary and can no longer influence the line for broader FS. 
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Fig. 10. The 525.0 nm line-core intensities (3-D NLTE left, LTE 
right) along the x-axis for three different FS atmospheres (solid 
HSRA (standard), dashed NET, dotted PLA). In the field free 
part, we always use HSRA. 



A few test computations with FS of different half-width, but 
with the same wall thickness also show a decrease in the core 
intensity with FT size. The decrease is only slightly less steep 
than in Fig. [9] confirming that the main effect is due to the FS 
size and not due to the wall thickness, as already suggested by a 
comparison of Fig. [9] with Fig. [8] 

3.6. Variation in the temperature stratification 

We then modified the model by replacing the at mospher e in the 
FS interior by, in turns, the network model of Solanki (1986) 
(from now on referred to as NET) and a plage model accord- 
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ing to lSolanki & Briglievkl (ll992h (hereafter PLA). The nominal 
field strength of the FS (B = 1600 G) and the external QS atmo- 
sphere (HSRA) were unchanged. Although these semi-empirical 
models are relatively old, th ey agree well w ith more recent 3-D 
radiation MHD simulations dVitas et alj|2009l) . The NET model 
is the hottest among the three selected atmospheres, being hotter 
than both the HSRA and the PLA model at all r. The PLA model 
is similar to the HSRA at heights around t c = 1, but has a flatter 
temperature gradient and is therefore hotter than the HSRA at 
the formation height of the 525.0 nm line core, which was taken 
into consideration here (see Fig.fTOb. 

Far away from the FS, the rest intensities are similar for all 
models and computational methods, which is unsurprising as the 
canopy lies so high that the FS atmosphere is almost irrelevant. 
Close to or inside the FS, however, the differences are significant 
for both the different models and methods. 

We first considered the LTE case, which is easy to un- 
derstand because the different temperatures of the inner atmo- 
spheres are directly visible in the resulting rest intensities. The 
closer to the FS we are, the more we see the influence of the FS 
atmospheres. In accordance with the higher temperature in the 
PLA and NET models, the rest intensities increase towards the 
FS (right panel in Fig.fTOb. The highest intensity is reached near 
the FS boundary and for the hottest model (NET). The weak- 
ening of Fe i lines in hotter atmospheres is due to the enhanced 
ionization of iron. 

In the NLTE case, the rest intensities vary strongly as well, 
but only in the interior of the FS. For the NET model, the NLTE 
intensity lies even below the LTE value. This can be explained 
by the large temperature difference between the NET model and 
the HSRA at the geometric height of line core formation. The 
downward shifted NET model in the FS is still approximately 
200 K hotter than the external HSRA model. We note that this 
no longer applies at greater depth. In the sheet, the LTE line- 
core intensity is therefore even higher than the corresponding 
3D NLTE intensity because the flux sheet interior is "cooled" 
by horizontal RT towards the now "cold" walls. The usual line 
weakening by the horizontal transport of UV radiation, as seen, 
e.g., in the HSRA-HSRA reference case, is inverted and the net 
horizontal transport of radiation has the opposite sign. 

The surprisingly good agreement between the PLA-LTE and 
the PLA-3D-NLTE rest intensity leads to the conclusion that the 
PLA model-when calculated in LTE-compensates the most ef- 
fectively for 3D and NLTE effects. The other lines considered 
in this paper also agr ee quite well (not plotted). This suggests 
that the PLA model of Solanki & Brig lievicl(ll992l) also remains 
valid as a description of the atmosphere of magnetic elements in 
full 3-D NLTE, at least for the spectral lines considered here. 

3. 7. Influence of geometry: Flux tubes 

We present results for a flux tube (FT) geometry, i.e. for mag- 
netic elements with a round cross-section. Despite the lower res- 
olution of the atmospheric model (cf. Sect. 12.3. U . the tube may 
be represented in a fairly smooth way because the filling fac- 
tor of 0.2 leads to a considerably larger FT fraction when we 
look along only one dimension (e.g. x-axis). The hot walls are 
therefore resolved almost as well as in the FS case. All other 
model ingredients (atomic model, atmospheric components, ra- 
dius, magnetic field, filling factor, and smoothing) remained un- 
changed from our standard FS model. 

Fig.QT|shows maps of the core intensities of three lines com- 
puted in 3-D NLTE (top) and LTE (bottom). The grayscales for 
the 3-D NLTE and the corresponding LTE images are the same 




I 



-200 -100 100 200 -2. 
x[km] 



LTE 
630.1nm 




1 1 .... m -70Q [ i , 

-200 -100 100 200 -200 -100 

x[km] xlkml 




m -200m 

200 -200 -100 100 200 



I 



Fig. 11. Intensity maps of a flux tube model in the cores of the 
three lines indicated above each panel. The grayscale for the 3- 
D NLTE (upper panels) and LTE (lower panels) cases are the 
same within the pair of plots for a particular spectral line. The 
scale varies, however, from one spectral line to another between 
zero intensity and the maximum achieved in 3-D NLTE (upper 
panels), which is set to 100%. 
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Fig. 12. Intensity in three line cores (from left to right: 525.0 
nm, 630. 1 nm, 630.2 nm) along a slice passing through the cen- 
ter of the rotationally symmetric flux-tube (FT) model. Vertical 
lines indicate the FT boundaries (solid) and the wall thickness 
(dashed). Solid lines 3-D NLTE; dashed lines LTE. 



for a given spectral line to facilitate the comparison. Owing to 
the rather different maximum intensities reached by the different 
lines, the greyscales were normalized individually. The maxi- 
mum intensity per line — typically reached near the FT boundary 
in the 3-D NLTE case — was used as the reference value (100%). 
Fig.Q~2]displays the same intensities along a cut passing through 
the axis of the FT. We observed the same effects as for the stan- 
dard flux sheet, although the line core intensities in the tube it- 
self are higher than for the FS, which may be explained by the 
more intense irradiation from the now completely surrounding 
hot walls. Differences near the FT boundary originated from the 
different geometrical expansion of the FT in three dimensions. 

This confirms that the results obtained for the (2-D) FS ge- 
ometry apply qualitatively to the 3-D FT case as well. The hori- 
zontal transfer effects, i.e. those of 3-D RT, are actually larger in 
FTs than for a FS (see e.g. the line weakening of the 525.0 nm 
line at the centers of the FT and FS). 

3.8. Center-to-limb variation 

Fig- 02] depicts the synthetic images of a flux tube in the inten- 
sity (3D-NLTE calculation) of the 525.0 nm line core for three 
different viewing angles. In the following, y = arccosyu de- 
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Fig. 13. Intensity map of the 3-D FT in the core of the 525.0 nm line for three different viewing directions as indicated on top of the 
panels. See Sect. l3.8l for definitions of y and*. The maximum intensity of the second image, showing the highest intensity values, 
was used to normalize the grayscale of all images. 
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Fig. 14. Core intensity of the 525.0 nm line for the different 
viewing angles considered in the 2-D FS geometry. Viewing an- 
gles and /i = cosy are given in the text labels. In a FS, rays with 
the same // (dotted, dashed) do not necessarily have similar in- 
tensities. The hot walls are more clearly seen for larger azimuthal 
angles *. Note that the location of the maximum intensity in x 
is — for inclined rays — somewhat arbitrary as the outgoing inten- 
sity is plotted at the top of the atmosphere. 



notes the angle between the vertical (z-axis) and the line of sight 
(LOS) and* the azimuthal angle in the xy plane measured from 
the x-axis. In the leftmost panel, showing the vertical view (i.e. 
[i - I .0), we see the strong and very localized (« 10-20 km) in- 
tensity maximum near the FT boundary and, less strikingly, but 
still clearly visible, the enhanced intensity in the interior of the 
FT. For ju = 0.93 (i.e. y = 21.5°) and* = 45° (middle panel), 
we still see half of the ring, which is now much broader (« 50 
km) and brighter. It actually resembles an arched structure such 
as a small loop. The rest of the hot wall is hidden from view by 
the optically thick cooler gas of the external atmosphere in the 
foreground. For// = 0.7 (y = 45°) and* = 0° (right panel), the 
hottest parts of the wall cannot be seen directly anymore, leading 
to a more diffuse intensity distribution, partially produced by the 
warm walls at greater height, partially due to the horizontal RT. 
We note that for all LTE calculations, no significant structures 



could be recognized at all. The same is true for NLTE computa- 
tions at n< 0.7 (y > 45°). 

Fig. [14] shows the intensity distributions for the same (and 
more) viewing angles, but now for the 2-D flux sheet. In the FS 
geometry, there may be different intensity distributions for a sin- 
gle fi and different azimuthal angles, as illustrated for// = 0.68. 
When the ray is mainly in the x-direction, i.e. * is small (dotted 
curve), we see the same smooth distribution as for the 3-D FT 
(rightmost panel of Fig.fTSl. If we look at the same fi = 0.68, but 
now mainly along the y-direction, i.e. x is large (dashed curve), 
we have a more direct view of one of the hot walls because the 
ray is oriented along the central plane of the flux sheet. Figure 
[14]shows an additional viewing angle with a very small // of 0.26 
(y = 75°, thin solid line), where one barely recognizes any vari- 
ation, since the hot wall is almost completely hidden by the front 
part of the FS. 

3.9. Solution for the full Stokes vector 

Here we discuss solutions for the full Stokes vector including 
Zeeman splitting (following the method described in Sect. 12.31 
and the Appendix in more detail). 

Figure Q3] displays the run of the line-core intensities of our 
four selected lines for different calculation methods. Each panel 
pair compares the results obtained by including Zeeman splitting 
(panels on the right) with those without it (panels on the left). We 
note that the non-magnetic case (Stokes / only) was already dis- 
cussed in Sect. 13.31 We now show both the 524.7 nm and 525.0 
nm lines as their profiles are no longer identical in the Zeeman 
split case due to their different Lande factors. 

The inclusion of Zeeman splitting may lead to either an in- 
tensification or attenuation of the NLTE and horizontal RT ef- 
fects, but the concrete behavior of a line can be understood in 
terms of its strength, Lande factor, and to a certain extent its 
temperature sensitivity. 

We first consider the 524.7 nm and 525.0 nm lines. In the 
absence of Zeeman splitting these two lines are practically iden- 
tical, as expected (panels a and c). In the panels b and d, however, 
the 525.0 nm line is more strongly weakened owing to its higher 
Lande factor. This is true for all three curves. The Zeeman split- 
ting essentially reduces the opacity in the line core. This may 
allow the hot walls to become visible to the line core, as happens 
in LTE for the 525.0 nm and 630.2 nm lines, where we observe 
strong peaks near the boundary in the Zeeman-split case. 

At the FS center, the effects of 3-D RT (i.e. the relative dif- 
ference between the 1-D and 3-D intensities) may also be ei- 
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Fig. 15. Core intensities for the four selected lines and different calculation methods (solid 3-D NLTE, dashed LTE, and plus 
signs 1-D NLTE; the vertical lines are as in previous figures). Each panel pair shows the intensities of a single line for a formal 
solution with Stokes / only (i.e. no Zeeman splitting; panels a, c, e, g) and one considering the full Stokes vector ( panels b, d, f, 
h) as indicated above each panel. Note the x-scale in the left panels of each pair (Stokes / only) is inverted to provide a clearer 
comparison of the intensities in the FS interior. 
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Fig. 16. Stokes / profiles of the 525.0 nm line at five selected x-positions for the standard FS model including Zeeman splitting. 
Solid lines 3-D NLTE, dashed lines LTE, and plus signs 1-D NLTE. Refer to Fig. [4] for the profiles calculated without Zeeman 
splitting. 



ther enhanced or damped. The 630. 1 nm line, which is the most 
Zeeman insensitive of the four considered spectral lines, displays 
little dependancy on Zeeman splitting in LTE and 1-D NLTE, 
but a strengthening of the 3-D effects. In contrast, precisely the 
opposite behavior is exhibited by the 630.2 nm line, where the 
relative difference between the three methods is very small when 
Zeeman splitting is included. The 524.7 nm and 525.0 nm lines 
are particularly affected by 3-D effects. We note that these curves 
refer only to the line core. At other places in the line profile, 
even the 630. 1 nm line appears to be considerably influenced by 
Zeeman splitting. The exact behavior depends on details of the 
splitting compared to the line width, as well as other factors. 

Figure [TBI demonstrates the variation in the Fei 525.0 nm 
Stokes / profiles with distance from the x-axis of the FS, in anal- 
ogy to Fig. [4] but for the solution of the full Stokes vector. As 
expected, far away from the FS, the profiles computed in 3-D 
are very similar to LTE profiles (panel e). The commonly made 
assumption that Fei 525.0 nm is "formed close to LTE" seems 
to be valid. The inclusion of the Zeeman effect has little influ- 
ence because the line is formed mainly in the non-magnetic at- 
mosphere below the canopy. The situation changes when we ap- 
proach the boundary of the FS where the field strength is higher 
and the line is partially formed in the canopy. In addition to the 
effects already discussed in Sect. 13.31 the inclusion of the full 
Stokes vector in the formal solution produces much broader pro- 
files and enhanced line-core intensities (especially in the 3-D 
NLTE and 1-D NLTE calculation, less so in LTE) at all locations 
close to or inside the FS (compare panels a) to c) of Fig. [Trjlwith 
the corresponding panels of Fig. |4j». 

Figure[17]displays the Stokes V/I c and Fig.[18]the Q/I c pro- 
files corresponding to the first four / profiles in Fig. [16] Positions 
further away from the FS do not show a significant polarization 
signals, as the line is formed in the field-free part of the atmo- 
sphere. 

The 3-D NLTE V/I c profiles are weaker than the LTE ones at 
all locations and, inside the FS, also weaker than the 1-D NLTE 
ones (except for the special case very close to the boundary, 
where the influence of horizontal RT is reversed). In the 525.0 
nm line, the larger Zeeman splitting in the 3-D profiles can be 
clearly seen (both in Stokes / and V). In general, however, ow- 
ing to 3-D RT, Stokes V displays smaller effects than Stokes / 
because these effects are largest in the line core, where Stokes 
V is zero. Interestingly, the form of the 525.0 nm V profiles in 
panel e) to differ from those calculated in either 1-D or LTE or 
those of the 524.7 nm line. The center of gravity of the 525.0 nm 
3-D profile shows a systematic shift away from the line center, 
an effect that may be of considerable relevance to the diagnos- 



Table 2. Wilson depression wq for different field strengths B 
at z — km (i.e. t c — 1 in the QS) and different FT model 
atmospheres (NET, PLA). 



B [Gauss] 


1700 


1600 


1500 


1400 


1300 


w D [km] 


258 


191 


149 


120 


97 


B [Gauss] 


1200 


1100 


1000 


NET: 1600 


PLA: 1600 


w D [km] 


79 


64 


51 


238 


206 



tics of the magnetic field strength, since such a profile mimics a 
stronger field strength and a seemingly lower filling factor than 
the profile in LTE. The Stokes V 3-D NLTE profile shape in the 
FS center is consistent with an enhanced weakening of the line 
core inside the FS. 

The Stokes Q/I c profiles plotted in Fig. [18] clearly illustrate 
that the effect of 3-D RT is largest near the line core. At the 
center of the FS, these profiles vanish for reasons of symmetry 
(the plot does not exactly corresponds to the FS center, so that 
the signals do not completely disappear). Near the flux sheet's 
boundary, however, as shown in panel b, the strong weakening 
of the 3-D NLTE Stokes Q in the line core is clearly apparent. In 
the line wings, unsurprisingly, the 1-D and 3-D profiles almost 
coincide. The difference in the strength of the 7r-components of 
the 524.6/525.0 nm lines is due to the Zeeman desaturation of 
the line core, which is stronger for the higher g line. 

3.10. Variation in the field strength 

The field strength has a two-fold effect. Firstly, it affects the con- 
struction of the atmospheric models by determining the Wilson 
depression, which in turn influences the magnitude of the hori- 
zontal transfer effects. The Wilson depression corresponding to 
the different field strengths is listed in Tab. [2] Obviously, the de- 
pendence is strongly non-linear. Secondly, the field strength di- 
rectly affects the line profile through the Zeeman effect. Here we 
considered the effect of changing the field strength, while other 
atmospheric parameters were unchanged. 

Figure [T9b illustrates the influence of the magnetic field 
strength on the 525.0 nm line core intensity without Zeeman 
splitting. As expected, the influence of the field strength is strong 
even when neglecting the Zeeman splitting. The main effects — 
the strong peak near the flux sheet's boundary and the enhanced 
intensity in the center of the FS — both increase with increasing 
Wilson depression. Below approximately 1500 G at z — km 
(which corresponds to t c — 1 in the HSRA), the peak due to the 
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Fig. 18. Same as FigQUbut for Stokes Q/I c . 



hot walls disappears, while the intensity produced by 3-D RT at 
the FS center remains somewhat higher than in 1 -D NLTE even 
for 1000 G, demonstrating that horizontal radiative transport is 
still a relevant factor. 

Figure [T9b depicts the intensities obtained when the line 
splitting due to the Zeeman effect is taken into account. The 



full Stokes solution provides significantly higher rest intensities 
close to and inside the FS compared to the solution neglecting 
the Zeeman effect, whereas the intensities far away are unaf- 
fected as they are mainly formed in the — non-magnetic — QS at- 
mosphere. We note that only part of the increase in rest intensity 
in Fig.[T9b is directly due to Zeeman splitting. Its indirect influ- 
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Fig. 19. Influence of the magnetic field strength on the 525.0 nm line-core intensity. Three-dimensional NLTE (solid, right part 
of each panel), 1-D (solid, left part of each panel), and LTE (dashed) for five different field strengths of 1700G (highest curves, 
thicker), 1600G (standard), 1500G, 1400G, and 1000G (lowest curves, thinner) at z — km. Panel a: Intensities achieved when 
Zeeman splitting is neglected. Panel b: those obtained after including the Zeeman effect. Vertical lines indicate the FS boundary 
(solid) and FS center (dashed). 



ence in terms of 3-D RT is equally important, as can be seen by 
comparing Fig.[T9h with[T9b. 

4. Discussion 

After presenting, in Sect. [3] the influence of various model pa- 
rameters on both the line profiles and the rest intensity, we now 
consider questions such as, how strongly do typical observa- 
tional quantities di ffer if one employs true 3-D RT instead of 1- 
D or even LTE? As IStenholm & Stenflol d 19771) already pointed 
out, not only the inclusion of NLTE effects, but also horizontal 
RT may play major roles when considering a flux tube geometry. 

Here we have gone beyond the earli er stud ies of 
Sten holm & Stenflol Jl977l)jCaccin & Severinol d!979l) . and 
Brul s & von der Ltihd d200lh by carrying out a systematic study 
of the changing parameters of a thin flux tube (or thin flux sheet) 
on the profiles of four widely used spectral lines at different 
levels of simplification in the RT (full 3-D NLTE, 1-D NLTE, 
and LTE). 

To compare the influence of different calculation methods 
and models on the determination of important atmospheric pa- 
rameters, such as temperature, abundances, magnetic flux, and 
field strength, we selected frequently used observables to quan- 
tify the effect of the different parameters. We show the results 
for the Fei 525. nm line, which is a line widely used in solar 
diagnostics (e.g., Bart hol et alJl2011h . The effects are similar for 
the other lines, though they may differ in their details. Zeeman 
splitting and UV opacity enhancement were included in all com- 
putations presented in this section. 

4.1. Core line depth and determination of temperature 

Table [3] gives an overview over the line depth di = (I c - r)/I c , 
where r is the rest intensity at line center and I c the continuum 
intensity, calculated for some of our model FSs. We consider 
here the depth of the spatially averaged profile, i.e. averaged over 
the area over which the flux sheet expands until it merges with its 
neighboring sheets. For a filling factor of 0.2, which is assumed 
for all computations presented in this section, this corresponds to 



Table 3. Line depth in units of l c of the spatially averaged Fe i 
525.0 nm line profile for 3-D NLTE, LTE, and 1-D NLTE com- 
putations, as well as for several model FSs derived from our stan- 
dard model. 



Model 


3-D NLTE 


LTE 


1-D NLTE 


Quiet Sun (HSRA) 


0.83 


0.87 


0.83 


B = 1000 G 


0.81 


0.86 


0.82 


B = 1300 G 


0.76 


0.85 


0.78 


B = 1400 G 


0.71 


0.84 


0.74 


B = 1500 G 


0.65 


0.81 


0.69 


B = 1600 G, standard 


0.59 


0.77 


0.63 


B = 1700 G 


0.53 


0.68 


0.56 


PL A (B = 1600 G) 


0.61 


0.62 


0.57 


NET (B = 1600 G) 


0.59 


0.58 


0.53 



integrating out to five times the distance from the FS boundary 
(at a height corresponding to r = 1 in the QS). 

The results obtained for d\ of Fei 525.0 nm are summa- 
rized below and in Table [3] The standard model corresponds to 
B = 1600 G, R = 100 km, a filling factor of 0.2, a wall thickness 
of 20 km, the HSRA in both the QS and the FS interior, and an 
enhanced UV opacity. The line depth obtained in 3-D NLTE in 
the standard model is printed in boldface in Table [3] The values 
for the QS entry were taken from a flat 1-D HSRA model with- 
out any magnetic structure. They correspond well to the values 
reached far away from the FS. We note that because of Zeeman 
splitting, the line core depths do not necessarily correspond to 
the maximum line depth. However, in the averaged profiles, the 
core intensity is always very close to the minimum intensity. The 
main conclusions reached from these computations are: 

1. The average 1-D profile is always deeper than the corre- 
sponding 3-D profile except for the hotter PLA and NET 
atmospheres. 

2. The values of d\ obtained for 1-D and 3-D RT differ by ap- 
proximately 5-10%. The LTE and 3-D line depths may differ 
by up to 30%. 
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3. The atmosphere describing the interior of the FS is — as 
expected — strongly relevant. The LTE and especially the 1- 
D NLTE line in the hotter NET atmosphere is weaker than 
that of the true 3-D NLTE calculation. We note that the dif- 
ference in d\ between the FS atmospheres is much smaller 
when true 3-D NLTE is used. 

4. For models with weak magnetic fields (< 1500 G), the effects 
of 3-D RT are relatively small owing to the small Wilson de- 
pression and the small amount of radiation entering the FT 
from the walls. For strong fields, the effects increase dispro- 
portionately. 

5. The influence of the variation in the geometric quantities 
(such as the flux sheet radius or the wall thickness, which 
are both not given in Tab. on the depth of the average line 
profile is generally small. 



Table 4. Same as Table [3] but for the equivalent width given in 
% of the nominal value obtained for 3-D NLTE in our standard 
model (boldface font). 



Model 


3-D NLTE 


LTE 


1-D NLTE 


Quiet Sun (HSRA) 


118 


129 


118 


B = 1000 G 


119 


132 


120 


B = 1300 G 


117 


133 


119 


B = 1400 G 


114 


133 


117 


B = 1500 G 


110 


133 


115 


B = 1600 G, standard 


100 


132 


110 


B = 1700 G 


84 


127 


99 


PL A (B = 1600 G) 


100 


107 


93 


NET (B = 1600 G) 


101 


99 


86 



In the past, the thermal structure o f the magnetic elements 
has generally been derived in LTE (e.g.lStenflolll975UCha pman 
19771 [19791 1 Walton! [l987t ISolankj[l986l; ISolanki & Briglievid 
19921) . To quantify the error in the temperature of the FS when 
the proper 3-D NLTE formalism was not applied, we tried to find 
a 1-D atmosphere that reproduces the 3-D NLTE line depth ob- 
tained for the standard FS model. As the 3-D NLTE line depth is 
weaker, we searched for a hotter model than the HSRA. By re- 
placing the HSRA in the FS interior with an atmosphere whose 
whole temperature curve is shifted up by a constant amount rela- 
tive to the HSRA, we found that the difference between the depth 
of the LTE and the 3-D NLTE line profiles may lead to an over- 
estimate of the temperature of approximately 50 K for a filling 
factor of 0.2. 

Modern telescopes reach higher spatial resolutio n and can 
begin to resolve ever smaller magnetic elements (iLagg et al.l 
2010). We therefore also considered the error in the temper- 
ature obtained from observations that fully resolve the mag- 
netic elements. If we considered the profile emerging along the 
center of the FS, then this temperature error could reach val- 
ues as high as 300 - 400 K (for the standard model). When 
compared with the 1-D NLTE profiles, the difference was only 
slight ly smaller, confirming th at 3-D RT is relevant as predicted 
by IStenholm & Stenflol d 1978b . Similar temperature differences 
of a few 100K were also confirmed by inversions (kindly pro- 
vided to us by Michiel van Noort) of the calculated spectra. As 
in the vast majority of the cases employed in the past, the inver- 
sions were based on LTE forward calculations. 

Although this may be an alarming result for those who deter- 
mine temperatures by comparing observations with LTE calcu- 
lations, one has to consider that our model is rather simple and 
may not correspond to a realistic scenario. Thus, we consider it 
important to investigate this question further using realistic 3-D 
radiation MHD simulations, which we plan to do in an upcoming 
publication. 

Interestingly, the line depth in 3-D NLTE for flux sheets 
whose interiors are described by the PLA and NET atmospheres 
is, however, rather similar to the line depths obtained in LTE for 
the same models. Since these models were deduced by employ- 
ing the LTE approximation, this is a heartening result, since it 
suggests that these models may also have at least some validity 
for a 3-D diagnostic RT. This might help explain the similarity 
of the PLA model to the temperature stratific ation in flux sheet s 
formed in 3-D radiation-MHD simulations dVitas et al.l 120091) . 
albeit the RT in the simulations is restricted to LTE. 



4.2. Equivalent width and determination of abundances 

The equivalent width is a widely used observable for determin- 
ing the abundances or empirical oscillator strengths. Table |4] 
summarizes the equivalent widths of the spatially averaged pro- 
files for the same set of calculations and with the same averaging 
method as in Sect. 14.11 All the values were normalized relative 
to the 3-D NLTE results for the standard model, which is conse- 
quently given a value of 100% (given in boldface). 

The typical equivalent width of the LTE lines is much larger 
than that of 3-D NLTE calculations (32% in our standard model). 
The 1-D NLTE value is about 5 to 10% larger. The LTE equiv- 
alent width is not strongly influenced by any variation in the 
model parameters because the LTE lines are typically saturated 
and do not vary strongly across our FS models. Exceptions are 
the calculations with hotter FS/FT atmospheres (PLA, NET), 
where the differences between the 3-D NLTE and LTE calcu- 
lation are much smaller. This is not unexpected given the results 
of Sect. O 

Assuming that the true equivalent width is closer to that ob- 
tained in 3-D NLTE, then the larger equivalent width obtained 
using LTE or 1-D NLTE can be misinterpreted as a lower iron 
abundance in the magnetic features or a higher temperature than 
really present, possi bly combined with reduced micro turbulence. 
ISheminova & Solankil (Q999) determined elemental abundances 
in magnetic elements under the assumption of LTE and found 
similar values as in the quiet Sun. However, they employed em- 
pirical model atmospheres for the magnetic elements that were 
derived in LTE, which likely had an incorrect temperature struc- 
ture to compensate for the wrong line depths obtained in LTE. 
The use of these models also gives line equivalent widths close 
to the observed values for the correct abundances. This implies 
that a self-consistent LTE computations should more often lead 
to errors in the determined temperature than the abundances. 

By using MHD simulatio ns with magnetic flu x typical of 
moderate active-region plage, Fabbian et al.l d2010h found evi- 
dence of smaller equivalent widths for the iron lines in the visi- 
ble, in good agreement with many previous results (e.g.[Splanki 
119931) . This reduced equivalent width is of relevance to the de- 
termination of abundances. However, they employed only LTE, 
so that they basically see the effect of the different temperature 
stratifications in plage and quiet regions. We expect that the use 
of 3-D NLTE could give rise to a sizeable additional correction. 

4.3. Contrast and feature size 

One of the interesting results of our calculations is the strong 
intensity peak near the FS/FT boundary seen, e.g., in Fig. [3] 
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Table 5. Same as Table [3] but for the maximum contrast in the 
Fei 525.0 nm line core (see main text for definition). Since the 
contrast is, by definition, zero for the quiet Sun, it has been re- 
moved. 



Model 


3-D NLTE 


LTE 


1-DNLTE 


B = 1000 G 


0.18 


0.04 


0.17 


B = 1300 G 


1.31 


0.19 


1.61 


B = 1400 G 


2.63 


0.49 


3.20 


B = 1500 G 


4.58 


1.26 


5.37 


B = 1600 G, standard 


7.06 


3.04 


7.83 


B = 1700 G 


9.64 


6.80 


10.51 


PLA (5 = 1600 G) 


6.32 


7.07 


9.00 


NET (B = 1600 G) 


6.18 


7.79 


9.27 



We now investigate the maximum contrast, defined as (I max — 
Iqs)IIqs, where Iq$ is the intensity in the quiet Sun and I max is 
the maximum intensity reached along the x-axis (generally near 
the boundary of the FS). We note that we assume an infinite spa- 
tial resolution of the observations here (in contrast to Sects. 14.11 
andO. 

Although we did n ot find such strong contrasts as 
ICaccin & Severing] (1 19791) . our results still show considerable in- 
tensity variation within the model. Table [5] gives the line-core 
contrast for the same model calculations as given in Sects. 14.11 
and 14.21 The contrast may be either large (strong fields, es- 
pecially in NLTE) or small (weak fields, especially for LTE), 
depending on whether the hot walls are in direct view of the 
observer. Not only do the LTE and NLTE calculations differ 
strongly: the 3-D and 1-D NLTE calculations may differ con- 
siderably as well as can be seen from Tab. [5] 

Another interesting finding is the very narrow geometri- 
cal width (^ 20 km) of th ese structures as pointed out by 
Bruls & von der Liihel (1200 ll) . If such strong transitions from the 
FS exterior to its interior exist, a resolution of at least 5-10 km 
is needed to see the maximum contrast. However, we stress that 
this test is no proof that such small structures will actually be 
observable, since horizontal energy transfer via radiation tends 
to smooth out the horizontal temperature gradients. Thus, when 
replacing the HSRA model in the FS by the empirically derived 
NET or PLA models, the large and narrow peaks in the line core 
intensity at the edge of the FS are either strongly damped or 
even disappear. This question warrants additional investigation 
with the help of 3-D radiation MHD simulations. 

4.4. Stokes profiles and determination of magnetic field 
strengths 

Figures [17] and [18] in Sect. 13.91 reveal that all Stokes parame- 
ters are strongly influenced by the calculation method and the 
input parameters. In Fig. [17] both the strengths of the V/I c pro- 
files and the profile shapes differ significantly between the LTE 
and the 3-D NLTE case. Both quantities directly influence the 
determination of the magnetic field strengths. In the present sec- 
tion, we test this us ing the so-called center of gravity method 
dRees&Semellll979l) . The results are given in Table [6] for the 
same models as in Table [5] 

For the models with the same atmosphere inside and outside 
the FS (HSRA/HSRA), the field strength obtained from the 3- 
D NLTE computations is approximately 15% higher than that 
deduced from the Stokes V profiles computed in LTE, except 
for the B = 1700G case, where the relative difference is larger. 



Table 6. Same as Tab.[5] but for the field strengths B cog , in Gauss, 
10 km from the FS center, determined by the center of gravity 
method in the Fe i 525.0 nm line. 



Model 


3-D NLTE 


LTE 


1-D NLTE 


B = 1000 G 


733 


615 


664 


B = 1300 G 


1054 


899 


973 


B = 1400 G 


1192 


1024 


1110 


B = 1500 G 


1369 


1178 


1280 


B = 1600 G, standard 


1627 


1392 


1518 


B = 1700 G 


2204 


1755 


1926 


PLA (B = 1600 G) 


1325 


1226 


1416 


NET (B = 1600 G) 


1349 


1322 


1412 



About half of this difference is due to horizontal RT effects, 
while the other half comes from NLTE (see the difference in B cog 
between LTE and 1-D NLTE in Table©. We note that the Stokes 
V profiles (Fig. [TTT i in LTE and 1-D NLTE are quite similar but 
that the Stokes / profiles (Fig. [Toll differ strongly and lead to this 
sizeable difference in B cog . 

Part of the difference in B cog between LTE and 3-D NLTE 
may be attributed to the different height of formation of the spec- 
tral line computed under these different assumptions. 

The use of different FS atmospheres also has a strong influ- 
ence: the field strength deduced from the NET model in LTE 
reaches almost the value obtained from the 3-D NLTE calcu- 
lation. For the slightly cooler PLA atmosphere, however, the 
difference in field strength between 3-D NLTE and LTE is dis- 
tinctly higher, although it remains significantly smaller than that 
found in models with the HSRA/HSRA combination. This result 
further supports the conclusion that for these two empirical FT 
models, calculations in LTE provide a reasonable approximation 
for more realistic 3-D NLTE computations. 

5. Conclusions 

With an improved version of the 3-D RH code (lUiten broek 
2000), we have calculated widely used diagnostic iron lines in 
various flux-sheet (FS) and flux-tube (FT) models. Our adaption 
of the formal solution to mon otonic, pa rabolic Bezier interpola- 
tion of the source function (lAuerll2003l) is crucial to preventing 
considerable errors in the calculated intensities and to achiev- 
ing better convergence. We have developed an enhanced Bezier 
interpolation scheme for the formal solution of the full Stokes 
vector problem (for details see the Appendices). 

We started from a standard expanding flux-sheet model, 
which consisted of a HSRA atmosphere in both the non- 
magnetic QS and the magnetic FS, a radius R = 100 km, field 
strength B = 1600 G at external r c = 1, wall thickness w — 20 
km, and a filling factor a = 0.2. We studied the behavior of four 
Fei spectral lines (524.7 nm, 525.0 nm, 630.1 nm, 630.2 nm) 
under three different levels of approximation for computing the 
radiative transfer (full 3-D NLTE RT, 1.5-D NLTE and 1.5-D 
LTE), each without and with the in clusion of Z eeman splitting 
(following the field- free method of lReeslll969l) . We considered 
the dependence of these lines on FS/FT parameters, such as ra- 
dius, field strength or Wilson depression, wall thickness, temper- 
ature stratification, and UV opacity. 

We confirm that 3-D RT has a significant effect on spec- 
tral line formation, in qu alitative agreement with the findings of 
IStenholm & Steimol < fl977b . who used a simple two-level atom i n 
a cylindrical flux-tube model, and lBruls & von der Liihel (1200 lh . 
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who employed 2-D calculations in a more realistic expanding 
flux-sheet geometry. In the following, we list our main conclu- 
sions. 

The use of a m ulti-level atom seems to s lightly increase 
the effects found by Stenhol m & Stenflol (Il977l) . T he inclusion 
of mi ssing UV opacity in the way proposed by iBruls et al] 
( 1992), however, strongly reduces these effects. Taken together, 
we found that the influe nce of 3-D RT is much sm aller than in- 
dicated by the results of Sten holm & Stenflol dl977l) . 

The expanding flux sheet/tube geometry can lead to large 
contrasts in the line core, particularly at loc ations where the hot 
walls of the flux sheet/tube become visible dCaccin & Severinol 
Il979h . For almost all the considered models, the NLTE compu- 
tations produce larger contrasts than LTE, The contrast (and the 
strength of the 3-D effects) increases strongly with increasing 
field strength (owing to the increasing depth of the Wilson de- 
pression). 

A striking feature of the emerging line radiation is the 
strong peak seen at the FS/FT boundary for many of the mod- 
els and spectral lines, where the flux tube walls are in the di- 
rect view of the observer. These peaks (bright rings in FTs) can 
be a lot narrower than the generally quoted 50 km horizontal 
mean-free-path for plane-parallel atmospheres even when com- 
puted in 3-D NLTE. If such sharp boundaries exist in the so- 
lar photosphere, than they should remain visible even in the 
pr esence of horizontal RT. This result agrees with the findings 
of IBruls & von derLiihd (120011) and promises a new discovery 
space for telescopes reaching a spatial resolution as fine as 10 
km in the solar photosphere. We note, however, that one rea- 
son why the FT or FS walls were so prominently visible is that 
they represent a sharp horizontal transition from optically thin 
to optically thick gas. We suspect that the sharp structures that 
do not display this transition become more smeared out by the 
radiation. 

The resulting intensity of a spherically symmetric FT model, 
calculated in 3-D, shows qualitatively the same effects as those 
of the computationally cheaper FS model. However, the influ- 
ence of 3-D RT is greater in the FT than the FS owing to the 
larger wall area radiating into a particular point in the interior of 
the FT. 

Three-dimensional RT generally leads to a weakening of the 
line in the FS/FT interior. However, the inverse effect also oc- 
curs if the interior is hotter than the exterior (at equal geomet- 
ric height), a situation occurring at certain heights if e mpirical 
models, such as the NET or PLA fl ux-tube atmospheres dSolankil 
1986: ISolanki & Brigrievidl 19921) . are used for the FS/FT inte- 
rior. 

The behavior of different lines may vary strongly. The most 
important factors are the depth of the Wilson depression, which 
influences the importance of horizontal RT, and the formation 
height of the line, which determines whether the radiation from 
the hot walls reaches the absorbing atoms directly. 

The average difference between an LTE and a full 3-D NLTE 
RT is on the order of 10% in the line depth and up to 30% in 
equivalent width. Hence, if LTE is assumed when determining 
the temperature in the FS/FT interior, significant errors can re- 
sult in the FT temperature when 3-D effects are not taken into 
account. The higher the spatial resolution, the larger the errors. 
For spatially resolved magnetic elements, the errors can reach 
300-400 K. Abundances may also be underestimated in LTE or 
1-D NLTE, if correct temperatures are used for the analysis. 

Non-LTE and 3-D RT effects are strong for all Stokes param- 
eters. The difference between Stokes profiles computed under 
LTE, 1-D NLTE, or 3-D NLTE is of similar relative magnitude 



to that for Stokes /. The application of either LTE or 1-D NLTE 
is thus expected to lead to systematic errors in the determina- 
tion of the magnetic field strength o n the order of from 1 0% to 
20% if the center of gravity method dRees & Semeli n~979) is ap- 
plied. Since much of this error is introduced by the line being 
formed at different heights in LTE and NLTE, this difference in 
field strength is equivalent to an error in the height to which a 
given field strength is assigned. 

By good fortune the empirical PLA and NET models turn 
out to give rather similar results in LTE as in 3-D NLTE, so that 
the validity of results obtained with these models is broader than 
previously thought. 

In summary, the present investigation provides evidence that 
the interpretation of observations of magnetic flux concentra- 
tions may change when 3-D NLTE RT effects are taken into 
account. Our quantitative estimates can be considered a rough 
guide only, because of the simplicity of the model employed 
here. Hence, the natural next step of this line of investigation 
is to employ 3-D MHD simulations and investigate commonly 
used diagnostics in both LTE and 3-D NLTE. Such an investiga- 
tion is planned. 
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Appendix A: Formal solution: Influence of the 
source function interpolation method 

In the short characteris tics (SC) method ( Ol son & Kun asz 1987; 
iKunasz & Auerl H988I) . the formal solution is integrated only 
along the "short" path from the nearest up- and downwind 
points. The incoming intensities — as well as all other quantities 
used for the formal integration at the up- and downwind points — 
have usually to be interpolated from neighboring grid points. 
The problem of polynomial interpolation at the up- and down- 
wind points, which may lead to spurious extrema, w as widely 
discussed in the literature until Auer & Paletou(1994) found an 
elegant and efficient solution using monotonic interpolation. 

In contrast, one finds very few references to the prob- 
lems that arise when the source function is interpolated with a 
parabolic Lagrange polynomial along t he short characteri stic as 
proposed and established as standard bv lAuer et alJ (Qj94). With 
S(t) — So + c\t + C2T 2 , the integral over S along the short char- 
acteristics is then given by the analytical expression 

hocal = CqSq + C\W\ + C 2 W 2 , (A.l) 

where // oeo ; is the locally produced intensity without irradia- 
tion from outside. The w, may be — in the case of Lagrange 
interpolati on — expressed as sim ple functions of known quan- 
tities (see Auer & Paletou 1994). The first ind ication of these 
interpolation problems is given by lAuerl (12003). who proposes 
monotonic Bezier interpolation instead of simple Lagrange in- 
terpolation to prevent spurious extrema of S along the integration 
path. 

To demonstrate that the proper interpolation of S is as im- 
portant as the correct interpolation of up- and downwind quan- 
tities, we present the situation of an overshooting source func- 
tion, which we call "a critical case", taken from our fluxsheet 
atmosphere simulation (the standard model, see Sect. 12.31 for a 
description). Fig. lA.ll shows the values of S at the up- and down- 
wind points (S i, S2) and the integration point (S 0). The unique 
Lagrange parabola through the three base points is given by the 
solid line. One immediately sees that S clearly takes values be- 
low the base points So and Si, while it should decrease mono- 
tonic ally from S to S l. From 

n 

/ = I ie - ATi + j S(T)e- T dT, (A.2) 


it directly follows that the resulting intensity must be wrong (i.e. 
too small), because the area under S is far too small. We note 
that the factor <?~ T does not change the qualitative result as it is 
on the order of unity over the relevant integration interval. 

Bezier parabolas are defined by two base point values and a 
single control point,which is always located at the intersection 
of the two tangents of S at the base points. This control point 
therefore "controls" the run of S between the two base points. 
The unique Lagrange parabola through So, Si, and S 2 is recov- 
ered in the interval from S to S 1 when the control point (cq 
in Fig. IA.U lies on the corresponding tangents of the Lagrange 
parabola (fine dotted). One imposes monotonicity in the inter- 
val from S to S 1 by constraining the control point to values 
between the two base point values of S . By moving the con- 
trol point up to the lower of the two values (c ), one gets a new 
parabola (that no longer passes through S 2!) with a monotonic 
run in the integration interval (dashed curve). Therefore, the area 
under S is no longer artificially reduced. The difference in the lo- 
cally produced intensity for our sample, taken from an existing 
calculation, is considerable. 
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Fig. A.l. Interpolation methods: The Lagrange parabola (solid) 
does not give a good representation of the run of S along the 
short characteristic from So to S \. Parabolic monotonic Bezier 
interpolation (dashed) is achieved by restricting the control point 
Co (defined by the two tangents (dotted) of each parabola) to 
So <= c' <- Si. 



We therefor e used monot onic parabolic Bezier interpolation 
as proposed by lAuerl d2003) for all our calculations. We per- 
formed tests with a snapshot taken from a realistic 3-D radiation 
hydro dynamic simulation run (i.e. for B — 0) with the MURAM 
code dVogler et al.ll2005l). Additionally, we us ed a 3-D version 
of the FAL-C model of iFontenla et alj d 1993b obtained by ar- 
ranging many identical 1-D atmospheres into a 3-D cube. In 
both cases, we found that about 3% of all formal solutions suffer 
from S overshooting in the integration interval. While the result- 
ing intensity for the FAL-C model is almost the same as for the 
Lagrange interpolation (most of the critical cases in FAL-C dis- 
play only quite weak overshooting because of the smoothness of 
the atmosphere), the resulting intensities in the snapshot may dif- 
fer considerably owing to the strongly varying atmospheric pa- 
rameters along the integration path. Therefore, the source func- 
tion may exhibit strong overshooting. 

An increasing number of publications have recognized this 
proble m in recent years. For instance, lHauschildt & Baron] 
(2006) used linear i nterpolation (instead of mon otonic Bezier) 
for critical cases and lKoesterke et al.l d2008), and Hayek (2008) 
used cubic Bezier polynomials. 

In more extreme atmospheres than our fluxsheet model (as 
found e.g. in realistic MHD simulations, which will be dis- 
cussed in future publications), the resulting intensity may even 
be wrong by orders of magnitude with disastrous consequences 
not only for the final result, but also for the convergence of the 
RT computations themselves. The convergence is affected be- 
cause the (usually diagonal) approximate Lambda operator used 
for the improvement of the iteration is generally interpolated in 
the same manner as the source function. 

In our experience, we have found that the monotonic, 
parabolic Bezier interpolation, which is identical to Lagrange in- 
terpolation when S does not overshoot in the integration interval, 



has much better convergence properties than the pure Lagrange 
interpolation, which can lead to a complete divergence for non- 
smooth atmospheres, such as those along rays passing through 
the walls of a flux tube. 

However, imposing either monotonicity or a special treat- 
ment in critical cases may lead to convergence problems as well 
because a change in the interpolation method at a certain point in 
the cube may influence neighboring points in such a way, that in 
the next iteration the interpolation method is changed back again 
to the original method. This may lead to an infinite flip-flop sit- 
uation, so that convergence is never achieved. We encountered 
these flip-flop situations in a few of our unsmoothed flux-sheet 
models. They occurred at a convergence level of a few percent 
in the population numbers and were much stronger when linear 
interpolation was used in the critical cases (a change from linear 
to parabolic may have more influence on the area under S than a 
change from monotonic Bezier to normal Lagrange parabolas). 

The problem of strongly overshooting source functions may 
also arise in the 1-D schemes as they interpolate S in the integra- 
tion interval with a Lagrange polynomial as well. In addition, we 
therefore implemented and used monotonic Bezier interpolation 
for our 1-D calculations. 



Appendix B: Formal solution for the full Stokes 
vector problem 

The solution to the formal integration of the RT equation for 
the full Stokes vector w ithin the short characteristics method of 
lOlson & Kunaszl d 1987b and iKunasz & Auerl d!988l) was given 
bv iRees et al.l (1 1989b . Their second method, which is based on 
a diagonal-element lambda operator (DELO) uses linear inter- 
polation of_Jhe_soyrrce_fiincdon_components in the integration 
interval. ISocas-Navarro et all d2000l) improved the accuracy of 
the source function integration by introducing parabolic interpo- 
lation. The RH code is based on this improved method. 

As mentioned in the previous section, the interpolation with 
a Lagrange parabola may cause great difficulty in the formal so- 
lution. However, these problems are not limited only to Stokes 
/. The overshooting of the parabola may affect the other Stokes 
components in the same way. We therefore used the monotonic 
Bezier interpolation for all components. 

Unfortunately, the full Stokes vector problem is even more 
intricate: the integration with a monotonic, Bezier-interpolated 
source function may still lead to unphysical results in special 
(but not uncommon) cases. Fig. lB.ll illustrates a pathological sit- 
uation where the integration of the source functions leads to the 
unphysical result of V > I. Owing to the different (in this case 
even opposite) curvatures of S / and S y (the latter being a rep- 
resentative of all Stokes components), S y is larger than S / for a 
large part of the integration interval (to to ti). The likelihood of 
this occurring is highest when S y is almost as large as S / at one 
point and much smaller (or negative) at neighboring point(s). 

We note that the problem arises because one always has to 
consider a certain neighborhood of the integration interval when 
determining the interpolating function (except for linear interpo- 
lation where the two end points are sufficient). This very general 
result of interpolation theory is also responsible for the over- 
shooting of S i in the case of parabolic Lagrange interpolation. If 
the neighborhood does not behave "well", it may produce unre- 
alistic behavior (i.e. oscillations or overshooting) in the integra- 
tion interval. 

The Bezier parabola helps to reduce the influence of these 
difficult neighborhoods as it is defined in a more local way, i.e. 
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Fig. B.l. Interpolation method for the full Stokes vector prob- 
lem. Even with monotonic Bezier interpolation, the source func- 
tion of a Stokes component, e.g. Sy (dashed thin) may become 
larger than S / (solid thick) in the range from So to S i . The stan- 
dard control point lies far above cP (thin plus signs, both 
defined by two tangents given by thin dotted lines) owing to the 
different curvatures of the source function parabolas. To prevent 
these situations, we chose the control point for the Stokes V com- 
ponent (Cq , thick plus sign), according to Cg . Our cl is de- 
fined by the ratio (S Jj v) - c ( ^)/(S f - cf) (with x equal to / or V), 
which has to be constant. This choice prevents the source func- 
tion integrals of the Stokes components from becoming larger 
than those of Stokes /. 



We note that in the DELO method of iRees et all (Il989l) . the 
integration of the source functions is only part of the solution. 
The determination of the full Stokes vector intensities requires 
the solution of a system of linear equations. Even with our in- 
tegration method, in very rare cases, the solution of the equa- 
tion system produces unphysical results with one or more Stokes 
components becoming larger than Stokes /. Consequently, we 
changed the RH code in such a way, that the polarized Stokes 
component fulfills the condition 

Ve 2 + U 2 + V 2 < I (B.2) 

after each formal solution. In all of these cases in which the re- 
lation was not fulfilled, the polarized components were reduced 
by a single factor / (i.e. Q new = fQ, U new = fU, V„ eK = fV) 
such that the equality in the equation above held, i.e. 

yOneJ 2 + Unen 1 + V neH 2 = I . (B.3) 



by the two end points and a local control point, thus limiting the 
influence of points outside the integration interval. We note that 
the control point is influenced by the neighborhood (i.e. S 2), but 
in a less strict way. By imposing monotonicity (i.e. the control 
point has to be between the two values of the endpoints), the 
influence of the neighborhood is reduced exactly for these cases 
where overshooting would otherwise take place. 

As all Stokes components contribute to Stokes /, it is rea- 
sonable to take the shape of S /(t) as a reference. If we replace 
c Q v) by cx}^ in such a way that the shape of S /(r) predefines 
the shape of S y(r) along the integration interval, the unphysical 
condition of S y > S j may be prevented. For Bezier parabolas, 
this can easily be achieved with 



MQUV) _ (QUV) (/) (J) 

J c i),new _ J _ C 

JQUV) c (QUV) ~ 17) T(7) 

C 0,new 1 C 1 



(B.l) 



In this way, we use only local information, i.e. the shape of S /(t) 
and the endpoint values of S quv, to define the run of S quv(t). 
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